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Using models of the quasi-chemical theory of solutions, the activity coefficients of silicon 
are calculated in the melts Fe-Si, Ni-Si, and Fe-Ni-Si. The calculated free energies of 
solution of liquid nickel and silicon in liquid iron in the interval 0 to 1400 kbar and 
1500 to 4000 K, shows that Fe-Ni-Si alloy is stable under the conditions of the outer 
core of the Earth and the cores of the terrestrial planets. The oxidation-reduction condi- 
tions are studied, and the fugacity of oxygen in the mantles of the planets and at the 
core-mantle boundary are calculated. The mechanism of reduction of silicon is analyzed 
over a broad interval of P and T. The interaction between the matter of the core and 
mantle is studied, resulting in the extraction of silicon from the mantle and its solution 
in the material of the core. It is concluded that silicon can enter into the composition of 
the outer core of the Earth and Venus, but probably does not enter into the composition 
of the cores of Mercury, Mars, and the Moon, if in fact the latter possesses one. 


The question of the composition and mech- 
anism of formation of the cores of the Earth 
and the terrestrial planets is still open. 
Among the various approaches to the solu- 
tion of this problem, one of the most signifi- 
cant is analysis of a number of processes 
from the physical-chemical standpoint. In the 
first approximation they can be reduced to 
the following questions: (1) what were the 
redox conditions in the protoplanetary cloud ; 
(2) what are the redox conditions in the 
mantle and at the core-mantle boundary; (3) 
what is the possibility of formation of alloys 
of iron with other chemical elements upon 
condensation from the protoplanetary cloud ; 
(4) what is the possibility of solution of 
chemical elements in the liquid iron core; 
and (5) what is the possibility that equilib- 
rium exists between the core and mantle. 

The experimental data currently available 
allow sufficiently accurate calculation of pro- 
cesses occurring under the P-T conditions 


corresponding to the interior of the Earth. 

After experimental studies (refs. 1 and 2) 
on the dynamic compressibility of iron-nickel 
alloy confirmed the hypothesis of Birch (ref. 
3) about the inconsistency of the concepts of 
formation of a core of this composition, it 
became conclusively clear that in addition to 
the basic iron-nickel components the core 
contains other, lighter substances which re- 
duce its density by 10 to 20 percent in com- 
parison to the density of nickel-containing 
iron under the same conditions. 

MacDonald and Knopoff (ref. 4) and Ring- 
wood (ref. 5) suggested that such substance 
might be silicon. 

It is the purpose of the present work to 
analyze the following problems: (1) the pos- 
sibility of solution of silicon and nickel in 
liquid iron over a broad interval of P and T; 
(2) the mechanism for inclusion of silicon 
in the cores of the planets due to its reduc- 
tion in the mantle; and (3) the mechanism 
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for inclusion of silicon in the cores of the 
planets due to interaction between the ma- 
terial of the core and mantle. 

Solubility of Silicon and Nickel 
in Liquid Iron 

Information on the inclusion of impurities 
in the core can be produced by calculating 
the free energy of solution of a chemical ele- 
ment or compound in liquid iron, which con- 
sists of the free energy of the process at a 
constant temperature and the contribution 
of high pressure: 

aGt ^ Gt + f AVdP 
J 0 

The first term in the right-hand portion of 
the equation can be obtained experimentally 
or found theoretically using well-founded 
models of the theory of solutions. The second 
term is the difference between the partial 
molar volume of a component in a solution 
and the volume of the pure component at the 
same pressure. 

THERMODYNAMIC PROPERTIES OF 
MELTS IN THE FE-NI-SI SYSTEM 
AT 1 ATM. 

Metals similar in their properties to iron 
form nearly ideal solutions (quasi-ideal solu- 


tions) with iron. 

Consider a dilute solution 

Ni(,) = [Ni] in Fe(„ 

For one mole of dissolved material, the 
change in partial free energy between its 
initial state in the form of a pure component 
and its final state in the form of an ideal solu- 
tion is 

aG't soln. — Gi — G°f — RTlncii = RTlnXi 

(3) 

where Xi is the molar fraction of the dis- 
solved material. 

The calculated values of aG^ (soln.) with a 
content of 10 wt. % Ni in the solution are 
presented in table 1. 

Liquid alloys of iron and silicon do not 
follow the regularities of ideal solutions, but 
are described in terms of stricter theories in 
which the close order, as well as the tempera- 
ture and concentration dependence of the en- 
ergies of the interatomic interaction, must 
be considered in melts (refs. 6 and 7). 

The activity of silicon in melts at various 
temperatures has been calculated by an asjm- 
metrical version of the quasi-chemical theory 
of solutions (ref. 7). The method of calcula- 
tion and initial thermodynamic character- 
istics were described in detail in reference 8. 

The activity of silicon in Fe-Si and Ni-Si 
melts at Xsi = 0.2, which corresponds to 11.2 
wt. % Si, was calculated by two methods. In 


lnffi 2 = InXs -h 1/RTX^, (XiAGr + X^AGf^) (1 + r, + X 2 -^) 

oJ(.2 


+ XI X^d + r,)l/RT(AGf^' - AGg'«) + V^z 


X 2 - XiT} + XtXz (1 + V 


+ 2X i(l 7,)ln(l + r,)+2X 1 X 2 ; 

_ 1 + T? ^ v (-Y 2 — .Yi) , 2(aG|*^> — aG? 
dXi 8X2 1 - yj X1X2 ZRT 

xrv-^i:^2(i + >7)2l!. 


(5) 


( 6 ) 
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the first case, calculations were made using 
the equation 


log rsi = 


7268 

T 


+ 0.71 


( 11 ) 


lna2 


.X^iQi2(1 + ■q) 
RT{l-rj) 


(Xi - X2,) + InXa + i/2Zlnr ~ -^17 

L ^2 

+ - X^r,) (4^) 

, XiXad + nV 

(Zi - X 2 ,) (X 2 - x^r,y 



(4) 


where t] is the close order parameter; Q '12 = 
Qi 2 + iZRT\n2) J2 takes into account the de- 
pendence of mixing energy on temperature, 
where 


- Q,2 = Air I -f- qXi = - aH\_ - qX^, 
q = aH^ — aH " ; Xi and X 2 

are the molar fractions of the first and sec- 
ond components. In the second case, equa- 
tions 5 and 6, found on the preceding page, 
were used. 

The value of the close order parameter rj, 
considering the deviation in distribution of 
atoms in a binary solution from random dis- 
tribution, can be determined from the equa- 
tion 


Let us now estimate the influence of 10 
percent Ni on ysi, o-si, and the free energy of 
solution of silicon in the ternary system Fe- 
Ni-Si. In the first approximation, this calcu- 
lation can be expressed by the equation 


log ysi Xpe 

(Pe-Ni-Sl) 


log ysi + .Xni 

(Sl-Fe) 


log ysi 

(Sl-Nl) 


( 12 ) 


Combining the earlier calculated data on 
ysi in liquid binary systems, we produce : 

logysi = - 4856 + 0.661 (13) 

(Sl-Pe-Ni) rp 

When 10 percent Ni is added to a melt, the 
values of ysi and Usi decrease. However, this 
reduction is slight and the differences in free 
energies of solution of liquid silicon in liquid 
iron and in liquid nickel-iron are also not 
great, amounting to approximately 1.2 kcal 
through the entire temperature interval. 

Calculation of the errors of thermodynamic 
parameters of Fe-Si, Ni-Si, and Fe-Ni-Si 
melts was described in an earlier work (ref. 
8 ). 


(B-l) (1 + 7 , 12 ) ^XiX2 - 57,12 = 0 (7) 

The probability parameter B was calculated 
by two methods : 

and 


„ r2(XiAGr + ^2 AGD1 

B = exp 

I ZRT 


where aG<®> = aG — aG<‘^'“'> , 


(9) 


and Z is the mean coordination number. 

Calculation of the averaged values of log 
ysi using equations (4) and (5) by the 
method of least squares in the interval 1500 
to 4000 K leads to the following equations. 

In the Fe-Si system 

log ysi= -^^+ 0.656 (10) 


INFLUENCE OF HIGH PRESSURE ON 
SOLUBILITY OF NICKEL AND SILL 
CON IN LIQUID IRON 

In order to calculate the free energy of 
the process of solution at high pressures, in- 
formation is required on the volumes of the 
pure liquid metals and the partial molar vol- 
umes of the metals dissolved in liquid iron. 

In earlier works (refs. 1, 2, and 9 through 
12) , results are presented from studies of the 
dynamic compressibility of pure metals, iron- 
nickel, and iron-silicon alloys up to megabar 
pressures. From these data, the partial molar 
volumes of Ni and Si were calculated at 
pressures of 300, 500, 1000, and 1400 kbar 
(ref. 8). 

The values of free energy of solution of 
liquid nickel (10 wt. %) and silicon (-.' 10 
wt. % ) in liquid iron as functions of P and T 
are presented in table 1. 


In the Ni-Si system 
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The free energy of solution of nickel de- 
creases with increasing temperature (this 
follows from equation (1) ) , since the partial 
molar volume of nickel is greater than the 
volume of the pure element, i.e., hV > 0. As 
composition changes, aG^ changes in a rather 
complex manner. 

Figure 1 shows the change in free energy 
of solution of liquid nickel in liquid iron as 
a function of concentration of Ni at various 
temperatures at a constant pressure of 1400 
kbar. The extremes on the curves result from 
the dependence of the partial molar volume 
of nickel in the melt on concentration of 
nickel. At parameters characteristic for the 
core-mantle boundary (-^ 1400 kbar), the 
free energy has a minimum, indicating the 
stability of the Fe-Ni melt with nickel con- 
centration 10-18%. 

The data produced indicate that up to 
pressures on the order of 500 kbar, nickel is 
dissolved in liquid iron at all of the tempera- 
tures studied. At pressures of 1000 to 1400 
kbar, its solution is possible if the tempera- 
ture is over 2000 K (table 1). 

Table 2 presents certain data that allow us 
to judge the distribution of temperatures and 
pressures in the terrestrial planets and the 
Moon. At the core-mantle boundary of Mer- 
cury, Mars, and the Moon\ the pressures do 
not exceed 300 kbar. The positive value of 
the contribution to the value of aGJ 

is not great; the free energy of solution is a 
negative value and nickel may therefore be 
included in the composition of the cores of 



Figure 1. — Change in free energy of solution of 
liquid nickel in liquid iron as a function of concen- 
tration of nickel in Fe-Ni system at HOO kbar and 
various Temperatures. 


' It is not clear if the Moon has a core. The latest 
data presented at the Soviet-American Conference 
on the Cosmochemistry of the Moon and Planets did 
not clear up the question. Our calculations show that 
iron-nickel alloy is stable under the P-T conditions 
of the interior of the Moon, but they neither study 
nor answer the question of the existence of a core. 


Table 2 . — Assumed P-T Conditions at the Core-Mantle Boundary for Earth-Type Planets 


' Planet 

Mean Radius of 
Planet 
(km) 

Mean Radius of 
Core 
(km) 

Pressure at 
Core-Mantle 
Boundary 
(kbar) 

Temperature at 
Core-Mantle 
Boundary 
(K) 

Mercury 

2437 

1700-2000 


1400-2000 

Venus 

6054 

~3000 


- 4000 

Earth 

6371 

3471 

-1400 

-4000 

Mars 

3386 

500-1000 

250-300 

1000-1500 

Moon 

1737 

200-300 

-50 

1000-1500 
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Figure 2a . — Dependence of oxygen fugacity on tem- 
perature in the mantle of the Earth. P = 200 

kbars. Reactions; 

1. FeSiO,(s) = SiO,(st.s) + Fe(l) + %0. 

2. FeSu.sOtfsp.s.) = V 2 SiO, ( st.s) + Fe(l) 
+ 

3. FeO(s) = Fe(l) + V 2 O. 

4 . FeSiO,(s) = FeO(s) + Si(l) + O* 

5. FeSu.sOi(sp.s) = FeO(s) + VzSifl) + 

6a, FeSio.sOi(sp.s) + MgSi>.sOi(sp.s) = MgSiO, 
(s) + Fe(l) + %0. 

6b, '2,FeSia.sOt(sp,s) = FeSiOi(s) + Fe(l) + V 2 O, 

6c. FeSu.sOt(sp.s) + MgO(s) = MgSu.iOt (sp.s) 
+ Fe(l) + 

6d. FeSiO,(s) + MgO(s) = MgSiO,(s) + Fe(l) 

+ %Oi 

7. ViFe,0,(s) = %Fe(l) +V 2 O. 

8. SiO,(st.s) = Si(l) + Ot 

these celestial bodies. At the core-mantle 
boundary of the Earth and Venus, the pres- 
sure is approximately 1400 kbar, but inso- 
far as the temperature at the boundary is. 
greater than 2000 K, the total change of the 
value aG^ < 0 and Fe-Ni cores are not stable. 

Thus, calculations show that Fe-Ni melts 
are stable in the ranges of temperatures and 
pressures characteristic for the terrestrial 
planets and the Moon. 


The free energy of solution of liquid sili- 
con and liquid iron is practically independent 
of temperature (due to the increase in activ- 
ity), but decreases quite sharply with in- 
creasing pressure (due to high negative 
volume change) , i.e., an increase in pressure 
facilitates the process of dissolution, table 1. 

Calculations show that throughout the en- 
tire interval of P, T, and concentrations, 
the aG^ for the solution of silicon in liquid 
iron or nickel-iron has a large negative value. 
This indicates the stability of Fe-Si and Fe- 
Ni-Si melts under the conditions of the outer 
core of the Earth and the P-T conditions of 
the cores of terrestrial planets. 

The stability of Fe-Si melts at the temper- 
atures and pressures of the core of the Earth 
and the other planets is a necessary but not 
sufficient condition since we do not know 
even the principal means and mechanism of 
inclusion of silicon in the core. In order to 
understand this, we must study the redox 
reactions in the mantle of the Earth and the 
other planets. 



Figfure 2b . — Dependence of oxygen fugacity on tem- 
perature in the mantle of the Earth. P = 600 kbars. 
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Redox Reactions in the Interiors 
of Terrestrial Planets 

The fugacity of gases in the mantle in the 
early stage of evolution of a planet (after its 
formation from the protoplanetary cloud, but 
before the descent of iron into the core) can 
be calculated on the basis of equilibria of me- 
tallic iron with ferromagnesium silicates and 
oxides. We will consider the fugacity of 
oxygen (/o^) as an indicator of the redox 
situation (refs. 13, 14, and 15). 

For the reaction aA + hB = cC -\- dD + 
O 2 , the nature of the change of fugacity as a 
function of phase composition, P, and T and 
is determined by the expression: 

log/oj, = log/Sj -f alogftx + bloga^ - 

clogftc — dlogao 

P 

- I f AVsdP 

2.Z0SRT j 

O 

where fo^ is the equilibrium fugacity of oxy- 



Figure 2c . — Dependence of oxygen fugacity on tem- 
perature in the mantle of the Earth. P = 1000 
kbars. 



Figure 2d . — Dependence of oxygen fugacity on tem- 
perature in the mantle of the Earth. P = HOO 
kbars. 


gen for a reaction in which all components 
are in the standard state; AVg is the change 
in volume of the condensed phases ; and is 
the activity of the iron-containing component 
in the solid solution. In the calculations we as- 
sume Oi = 0.1. All calculations were per- 
formed with consideration of compressibility 
data. 

The I’esults of calculating the fugacity of 
oxygen for a number of limiting buffer re- 
actions in the interval 200 to 1400 kbar and 
2000 to 4000 K are presented in figures 2a, 
2b, 2c, and 2d. 

EARTH, VENUS 

Interpretation of the results of calcula- 
tions is based on available information on 
the composition and structure of the con- 
centric layers of the Earth, obtained from 
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experimental laboratory data at high pres- 
sures and seismology. Information on the 
structure of the other planets is quite sparse. 
Therefore, analogy is assumed between the 
chemical composition of the Earth and the 
terrestrial planets ; it is also assumed that all 
planets have iron cores. 

Let us study the limiting buffer reactions 
which controlled the fugacity of oxygen at 
various depths in the Earth during the early 
stages of its development. 

At 200 kbar 600 km), the lower limit 
of /02 was controlled by reaction (2) for 
reduction of spinel to Fe and stishovite (fig. 
2a) . At lower values of fo^, ferromagnesium 
silicates are reduced to Fe -|- Si02 or FeO + 
Si, contradictory to the seismic data which 
indicate that the transition layer of the man- 
tle consists of Fe and Mg silicates. 

The upper limit of /o, is probably con- 
trolled by the equilibrium of iron and mag- 
netite, while the optimal value of /o^ is that 
which changes according to curve (6) which 
reflects the Fe-Mg, silicate-Fe equilibrium 
and essentially combines four reactions (6a, 
b, c, d). 

We can see from figure 2a that the reduc- 
tion of liquid silicon and its formation of a 
melt with iron can occur under conditions 
corresponding to the intersection of reactions 
(2) and (5), i.e., at log > — 2.8 and T 
> 3200 K. 

However, these temperatures are signifi- 
cantly higher than the melting point of iron 
and of silicate matter (refs. 16, 17, and 18) 
and, judging from known thermal models, 
planetary interiors probably never reached 
such values. Consequently, the reduction of 
silicates to silicon at these depths could not 
have occurred. 

Calculations show that the process of re- 
duction becomes possible at pressures of over 
400 kbar (--' 1000 km) . 

The Birch-Magnitsky hypothesis presumes 
that the lower mantle basically consists of 
the oxides FeO and MgO (NaCl-type struc- 
ture) and SiOo (stishovite). Recently, this 
has been experimentally confirmed in experi- 
ments on the decomposition of Mg 2 Si 04 to 
periclase and stishovite at 330 kbar (ref. 19) . 


We can therefore accept the buffer equilib- 
rium between stishovite and silicon as the 
lower limit of in the lower mantle of 
the Earth, although we cannot exclude the 
possibility of conversion of the rutile-like 
modification of SiOo into the fluorite-like 
modification and we can accept the Fe-Fe 304 
equilibrium (7) as the upper limit (figs. 2b, 
2c, and 2d). 

The mass of the current terrestrial core is 
about one-third of the mass of the planet; 
consequently, it is obvious that the system 
Fe-FeO should have controlled in the 
lower mantle of the Earth and of any other 
planet until the iron sank into the core. In 
other words, this system was responsible for 
the distribution of oxygen in the mantle, and 
all other systems were of secondary signifi- 
cance. 

At a pressure of 600 kbar ( ^ 1400 km) , 
Tm (Fe)'=3100K (ref. 16). The reduction 
of stishovite becomes possible at T > 3500 K 
and log /©2 >3.1 (fig. 2b) . Under these con- 
ditions, liquid silicon is dissolved in the iron 
core forming an Fe-Si melt, which later 
forms the outer core. At 1000 kbar 2200 
km) and 1400 kbar (^ 2900 km), (Fe) 

3600 K and 4100 K, and the reduction of 
silicon could have occurred at T > 3900 K 
and 400 K (fig. 2c, 2d) . Calculations of the 
thermal history of the Earth show that such 
temperatures were reached during the first 
hundreds of millions of years of its evolution. 

Calculations of the thermal history of Ve- 
nus (refs. 20, 21, and 22) show that the dis- 
tribution of temperatures in the interior of 
Venus differs little from the terrestrial values 
and allows us to assume the presence of an 
iron core (ref. 23). We will probably not be 
far off if we assume that the reduction of 
silicon in the core of Venus occurs at the 
same P-T parameters as on Earth. 

As was pointed out in the first section, the 
free energy of solution of silicon in iron or 
nickel-containing iron is negative (table 1), 
which indicates the formation of Fe-Ni-Si 
melts in the lower mantle of the Earth and 
Venus. If we accept the scheme of Elsasser 
(ref. 24), this melt has an unstable equilib- 
rium and descends, forming the outer core. 
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Figure 3 . — Oxygen fugacity as a function of temper- 
ature at the presumed core-mantle boundanry of 
Mars. 


MARS 

Calculations of models of Mars with an 
iron core show that the core amounts to not 
over 4 percent (ref. 25) or 9 to 10 percent 
(ref. 26) of the mass of the planet. Calcula- 
tions of the thermal history of Mars, as a 
function of the concentration of radioactive 
elements, show that either the interior of the 
planet below 300 km is at the present time 
partially or fully molten or 50 percent of the 
volume of Mars is in the molten state, while 
the central parts remain solid (ref. 27). In 
the opinion of Sharma (ref. 28), modern 
Mars should have lower temperatures than 
the Earth at the same depths and should 
contain no liquid core. On the other hand, 
calculations of the thermal history performed 
by Toksoz (ref. 23) indicate a solid mantle 
and molten core. 

The pressure at the core-mantle boundary 
of Mars is approximately 300 kbar (refs. 29 
and 30). 


Figures 2a and 3 show the results of cal- 
culation of the fugacity of oxygen for a num- 
ber of reactions at 200 and 300 kbar. 

At 200 kbar, the mantle of Mars consists 
of solid solutions of silicates and, as was 
noted earlier, the lower limit of reduction of 
silicon corresponds to a temperature of 
> 3200 K. 

At a pressure of 300 kbar, two versions 
are possible: either a silicate or an oxide 
mantle. In the former case, the reduction of 
silicates and formation of Fe-Si melt is pos- 
sible at log /02 > — 1.2 and T > 3650 K (at 
the intersection of reactions 2 and 5) , the sec- 
ond case is possible where log > — 0.5 
and T > 3200 K (at the intersection of reac- 
tions 3 and 8). 

Study of the thermal history of Mars does 
not allow us to judge whether such tempera- 
tures were reached in the early stages of its 
development. Processing of data produced by 
Mariner 6 and 7 allowed Anderson (ref. 31) 
to propose that the entire core of the planet 
was not separated from the mantle. This 
means that the heating of Mars was insuffi- 
cient to melt all of the available iron, of 
which the core is made. Analogy with the 
Earth also gives us reason to believe that 
temperatures of 3200 to 3600 K have not 
been reached in the interior of Mars. 

Thus, using the currently available infor- 
mation on the thermal history and models 
of Mars, we can conclude that the reduc- 
tion of silicates and stishovite to silicon has 
not occurred in the interior of the planet and, 
consequently, silicon should not enter into 
the composition of the core. 

MERCURY, MOON 

The high mean density of Mercury indi- 
cates that it has an iron core amounting to 
some 60 to 70 percent of the mass of the 
planet (refs. 25 and 32). According to the 
calculations of Plagemann (ref. 33), the core 
of Mercury has a radius of 2112 km. Process- 
ing of the data from Mariner 10 has indi- 
cated that Mercury has a magnetic field and 
that the presumed diameter of the core is 
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rather large, approximately equal to that of 
the Moon (ref. 34). 

Iron is practically free of radioactive ele- 
ments, and, therefore, their concentration in 
the material of Mercury should be signifi- 
cantly less than in the material of other 
planets. Calculations of the thermal history 
with an initial temperature at the center of 
1000 K show that the temperature in the in- 
terior never rose above 2300 K, which is 
insufficient for melting and the formation of 
a core (ref. 20). On the other hand, with a 
higher initial temperature, say 1100 K, we 
find that the core might be either liquid or 
solid (ref. 23). 

In spite of the many studies of specimens 
of lunar soil, measurements of heat flux, 
and interpretation of seismic data, the exis- 
tence of a lunar core is still problematical. 
If a core exists, its radius does not exceed 
200 to 300 km, with a temperature in the 
core of about 1500° C (ref. 35) . 

The characteristics of the material compo- 
sition of the Moon and Mercury are not clear. 
We can assume, by analogy with Earth, that 
their interiors consist of ferromagnesium 
silicates. The pressure at the center of the 
Moon is about 46 kbar; at the core-mantle 
boundary of Mercury about 100 kbar (ref. 
25). As the available calculations show, the 
reduction of silicates to silicon can occur at 
such pressures at temperatures of over 2500 
to 3000 K. 

Since these temperatures are apparently 
unrealistic for these celestial bodies, we can 
conclude that silicon should not enter in the 
composition of their cores. 

Anderson (ref. 31), based on thermody- 
namic calculations of the condensation of the 
protoplanetary cloud, states his opinion that 
the Moon accumulated from Al, Ca, and Ti 
compounds such as gehlenite (CaAl 2 SiOT), 
spinel (MgAl 204 ), and perovskite (CaTiOs), 
which formed the primary condensate; it 
also cannot be excluded that Mercury is rich 
in refractory compounds, in addition to iron, 
and poor in the silicates of iron and mag- 
nesium. However, in this case silicon also 
should not be included in the composition of 
the core. 


In summary, as a result of processes of 
reduction and high temperatures in the lower 
mantles of Earth and Venus in the early 
stage of their evolution, silicon might have 
dissolved in the liquid iron or nickel-iron, 
and this alloy might form the outer cores of 
the planets. In the interior of Mars, Mercury, 
and the Moon, the temperatures were insuffi- 
cient for reduction of silicon, and it should 
not be included in the cores of these celestial 
bodies. 

In preceding sections, we have studied the 
possibility of formation of an iron-silicon 
melt, which forms the cores of planets in the 
initial period of their evolution. 

However, these processes do not limit the 
possible entry of silicon into the composition 
of the core. There is also another possibility: 
namely, the entry of silicon into the composi- 
tion of the core as a result of interaction of 
the matter of the core and mantle, which 
exists at the present time. 

Interaction Between Core and 
Mantle 

Seismic data show that at the present time 
the mantle of the Earth is in the solid state 
and its outer core is in the liquid state. The 
interaction between the matter of the mantle 
and core can be represented in the form of 
the interaction between stishovite and liquid 
iron or nickel-iron: 

Si02 (st. 5) + 2Fe(l) = 2FeO(5)4- [Si] 

where [Si] represents a silicon solution (11.2 
wt. % ) in liquid iron or nickel-iron. 

The values of free energies of solution of 
silicon are presented in table 1 ; the free en- 
ergies of reaction are calculated according to 
the data of Tret’yakov (ref. 36), Robie and 
Waldbaum (ref. 37), and Kubashewski et al., 
(ref. 38). 

The results of calculating the free energy 
of the reaction (AGp) in the interval 2000 
to 4000 K and 300 to 1400 kbar are presented 
in figure 4. 

Free energy is a negative quantity in the 
pressure interval in question if the tempera- 
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Figure 4 . — Change in free energy of the reaction 
SiO,(st.s) + 2Fe(l) = ZFeO(s) + [SiJCm Fe(l) 
as a function of pressure at 2000, 3000, and 
iOOO K. 

ture is over 3000 K, i.e., in this case the 
reaction will be shifted to the right. The 
pressure at the core-mantle boundary of the 
Earth and Venus is ^ 1400 kbar, the temper- 
ature ^ 4000 K, and the free energy of the 
reaction — 55.3 kcal/mol. Under these con- 
ditions, the core extracts silicon from the 
mantle ; this may explain its growth. 

The temperature in the interior of Mars 
is probably not over 1500 to 2000° C (ref. 
26), while the pressure at the core-mantle 
boundary is approximately 250 to 300 kbar. 
Therefore, if its mantle consists of oxides, 
the temperatures are insufficient for the oc- 
currence of this reaction. 

The pressure in the interior of the Moon 
and Mercury is insufficient for the formation 
of stishovite (but it is possible that coesite 


is formed) , and the mechanism of interaction 
between the core and the mantle just studied 
is unsuitable for these bodies. 


Conclusions 

There are therefore two possible versions 
of inclusion of silicon in the composition of 
the iron-nickel cores of the planets. The first 
proposes the reduction of stishovite, with 
subsequent formation of an iron-silicon melt 
which, due to the unstable equilibrium in the 
gravity field of the planet, sinks to the cen- 
ter. The second possibility is found in the 
current interaction of the core and mantle, 
as a result of which silicon is extracted from 
the mantle and dissolved in the material of 
the core. Both versions are realized only for 
the Earth and Venus. The temperatures and 
pressures in the interior of Mars, Mercury, 
and the Moon are insufficient for inclusion 
of silicon in their cores. It is possible that the 
inclusion of other elements must be studied in 
this case, for example sulfur or carbon. Only 
strict quantitative calculations can serve as 
proof. 

References 

1. Al’tshuler, L. V., G. V. Simakov and R. F. 

Trunin, Izv. AN SSSR, Earth Physics Series, 
No. 1, 1968. 

2. McQueen, R. G. and S. P. Marsh, J. Geophys. 

Res., Vol. 71, No. 6, 1966. 

3. Birch, F., J. Geophys. Res., Vol. 67, 1952, p, 227. 

4. MacDonald, G. J. F. and L. Knopoff, The 

Geophys. J. Roy. Astron. Soc. Geophys., Sup- 
plement M. N. R.A.S., Vol. 1, No. 4, 1958. 

5. Ringwood, a. E., Geochimica et Cosmochimica 

Acta, Vol. 15, No. 4, 1959. 

6. Vagner, K., The Thermodynamics of Alloys, 

Metallurgizdat Press, Moscow, 1957. 

7. Gel’d, L. V., B. A., Baum and M. S. Petrushev- 

SKIY, Ferroalloy Production Melts, Metallur- 
giya Press, Moscow, 1973. 

8. Kuskov, O. L., Geokhimiya, No. 11, 1974. 

9. Kormer, S. B. and a. I. Funtikov, Izv. AN 

SSSR, Earth Physics Series, No. 6, 1965. 

10. Pavlovskiy, M. N., Fiz. Tv. Tela., Vol. 9, No. 2, 

1967. 

11. Clark, S., Handbook of the Physical Constants 

of Rock, Mir. Press, Moscow, 1969. 



242 


COSMOCHEMISTRY OF THE MOON AND PLANETS 


12. Balchan, a. S., and G. R. Cowan, J. Geophys. 

Res., Vol. 71, No. 14, 1966. 

13. Khitarov, N. I. AND 0. L. Kuskov, lASPEl 

Symposium, Lima, Peru, 1973. 

14. Khitarov, N. I. and O. L. Kuskov, In : Solutions 

Solides en Mineralogie, Vol. 1, Colloque Inter- 
national, Orleans, Prance, May 27-30, 1974. 

15. Kuskov, O. L. and N. I. Khitarov, Geokhimiya, 

No. 1, 1974. 

16. Zharkov, V. N., Izv. AN SSSR, Geophysics 

Series, No. 3, 1959. 

17. Uffen, R. J., Trans. Am. Geophys. Union, Vol. 

33, No. 6, 1952. 

18. Higgins, G. and G. C. Kennedy, J. Geophys. 

Vol. 76, No. 8, 1971. 

19. Kumazawa, M. et al.. Nature, Vol. 247, No. 5440, 

1974. 

20. Mayeva, S. V., Physics of the Moon and Planets, 

Nauka Press, Moscow, 1972. 

21. Lee, W. H. K., Earth Planet. Sci. Letters, Vol. 4, 

No. 4, 1968. 

22. Pricker, P. E. and R. T. Reynolds, Icarus, Vol. 

9, No. 2, 1968. 

23. Toksoz, M. N., this collection. 

24. Elsasser, W. M., Early History of the Earth. 

In : Earth Science and Meteorites, 1963. 


25. Kozlovskaya, S. V., Physics of the Moon and 

Planets, Nauka Press, Moscow, 1972. 

26. Binder, A. B. and D. R. Davis, Phys. Earth 

Planet. Interiors, Vol. 7, No. 4, 1973. 

27. Mayeva, S. V., Dokl. AN SSSR, Vol. 159, No. 3, 

1964. 

28. Sharma, j. P., j. Geophys. Res., Vol. 78, No. 23, 

1973. 

29. Kozlovskaya, S. V., Astron. Zhum., Vol. 43, 

No. 5, 1966. 

30. Koslovskaya, S. V., Izv. AN SSSR, Earth 

Physics Series, No. 7, 1972. 

31. Anderson, D. L., The Moon, Vol. 8, Nos. 1 and 2, 

1973. 

32. Reynolds, R. T. and A. E. Summers, J. Geo- 

phys. Res., Vol. 74, No. 10, 1969. 

33. Plagemann, S., j. Geophys. Res., Vol. 70, No. 4, 

1965. 

34. Ness, N. R. et al., this collection. 

35. Lathem, G. et al., this collection. 

36. Tret’yakov, Yu. D., The Thermodynamics of 

Ferrites, Khimiya Press, Leningrad, 1967. 

37. Robie, R. a. and D. R. Waldbaum, U.S. Geol. 

Survey Bull., No. 1259, 1968. 

38. Kubashewski, O., E. L. Evans and C. B. Al- 

COCK, Metallurgical Thermochemistry, 1967. 



